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ABSTRACT 

A  procedure  for  the  design  of  a  recursive  digital  notch- 
filter  with  one  coefficient  is  proposed.   This  type  of  filter 
can  be  implemented  readily  with  CCD  (charge-coupled  device) . 
The  effect  of  coefficient  on  notch-gain,  notch-width,  and 
passband-ripple  is  examined.   The  computer  implementation  is 
tested  for  second,  third,  and  sixth  order  notch-filters  and 
the  optimum  coefficient  is  found  based  upon  limited  computer 
wordlength . 
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I.   INTRODUCTION 

Digital  filters  are  finding  wide  applications  in  many 
diverse  areas,  such  as  speech  processing,  picture  process- 
ing, sonar  and  radar  system,  digital  control  system,  etc. 
A  very  useful  type  of  such  filters  is  the  notch-filter 
which  attenuates  highly  a  particular  frequency  component 
in  the  input  signal  while  leaving  nearby  frequencies 
relatively  unattenuated. 

A  linear  digital  filter  can  be  realized  physically  in 
terms  of  the  three  basic  operations  of  addition,  multipli- 
cation and  unit  delay.   The  advantages  gained  by  realizing 
a  digital  filter  rather  than  its  continuous  counterpart 
lie  principally  in  its  stability,  accuracy,  reliability 
and  flexibility  [Refs.  7  and  20].   Furthermore,  its 
frequency  characteristics  can  be  varied  without  altering 
the  hardware  simply  by  changing  the  stored  set  of  filter 
coefficients.   The  disadvantages  of  digital  approach  lie 
in  two  extra  design  considerations.   Firstly,  the  choice 
of  wordlength  of  the  computer  and  analog- to-digital  con- 
verter.  These  problems  do  not  arise  in  the  corresponding 
continuous  design.   A  design  filter  operates  on  an  input- 
sampled  signal  by  means  of  a  computational  algorithm.   It 
can  be  simulated  on  a  general  purpose  computer  or  can  be 
constructed  with  special  purpose  digital  hardware,  such 
as  CCD  (charge-coupled  device)  implementation;  this  part 


is  beyond  the  the  scope  of  research.   Due  to  computer 
limited  wordlength,  the  round-off  error  in  computer  program 
is  produced.   This  can  cause  the  filter  characteristics  to 
deviate  considerably  from  the  ideal. 

This  research  attempts  to  present  a  survey  and 
discussion  of  the  effect  of  the  finite  wordlength  on  the 
accuracy  of  digital  notch-filter.   The  design  procedure  of 
digital  notch-filter  was  to  propose  rational  functions  of 
a  discrete  variable  which,  when  transformed  into  the  fre- 
quency domain,  have  desired  filtering  characteristics.   The 
design  was  proceeded  directly  in  terms  of  the  variable,  z, 
associated  with  discrete  time  intervals  rather  than 
through  the  intermediate  variables  associated  with  con- 
tinuous time. 

The  digital  notch-filter  with  one  coefficient  was 
investigated  and  tested  by  implementing  on  computer  IBM 
360/70  installed  at  the  Naval  Postgraduate  School. 


II .   THE  CHARACTERISTIC  TRANSFER  FUNCTION 
OF  DIGITAL  NOTCH-FILTER 

A.   TRANSFER  FUNCTION 

As  the  first  step  in  the  design  of  the  digital  notch- 
filter,  one  often  uses  the  bilinear  transformation  to 
approximate  the  analog  filters  in  discrete  time.   This 
transformation  process  is  perhaps  the  most  utilized  approach 
relating  the  continuous  and  discrete  domains.   The  approach 
in  this  paper  differs  by  starting  directly  in  the  discrete 
domain  without  the  bilinear  transformation. 

The  approach  is  to  uniformly  sample  the  signal  u(t)  every 
T  seconds,  and  use  the  resulting  sequence  as  the  input  to  a 
recursive  digital  filter  governed  by  the  general  difference 
equation 


v(K)  =  bn   u(K)  +  bn  u(K-l)  +  ...  +  b   u(K-m) 


ax  v(K-l)  -  a2  v(K-2)  -  ...  -  amv(K-m)     (2-1) 


where  u(k)  and  v(k)  denote  the  values  of  the  filter's  input 
and  output  signals,  respectively,  at  the  K-th  instant. 

For  minimum  noise  generating,  cost,  and  simplicity,  one 
can  reduce  equation  (2-1)  to  the  minimum  number  of  coeffi- 
cients such  that  (2-1)  is  still  in  the  recursive  form. 


v(K)  =  u(K)  +  b   u(K-m)  -  a   v(K-m)         (2-2) 


The  above  simple  recursive  equation  of  the  recursive 
digital  filter  is  derived  from  (2-1)  by  letting 


and 


bo  -  l 


b,  =  b_  =  . . .  =  b   ,  =  0 
12  m-1 


a,    a_    •••    a   «  -  u 

1    2  m-1 


Then  its  corresponding  transfer  function  has  the  form 


H(z)  = 


1  +  b  z 

m 


-m 


1  +  amz 
m 


m 


(2-3) 


For  a  notch-filter  one  wishes  the  filter  characterized  by 
(2-3)  to  effectively  transmit  all  sinusoidal  inputs  with  the 
exception  of  the  wn  radiant/sec.  sinusoid  which  it  should 
completely  reject.   Unfortunately,  it  is  impossible  to 
implement  an  ideal  notch-filter 


Figure  (2-1)   Gain  frequency  behavior  of 
an  ideal  notch  filter. 
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B.   GAIN  FACTOR  OF  A  DIGITAL  FILTER 

The  poles  and  zeros  of  the  transfer  function  (2-3)  are 
computed  as  follows 


zm  +  b   =  0 
m 


.  my-, —  _  m  /r—  i  (2n+l) 

hence    z  =   /-b   =   /b   eJ  v    J 

m      m 


and     z   +  a_  =  0 
m 


n=0,l,2, 


m  . m  . —   i  f2n+ll- 


m  r-   J  (2n+l)- 
/-a   =   /a  eJ  v    ym 
m      m 


The  zeros  and  poles  are  respectively  given  by 


in 


m 


(2-4) 


and 


m, —  j(2n+l)- 

p.  =   /a_  eJ      ^m 


m 


i  =  l , .  .  .  ,m 


(2-5) 


From  (2-4)  and  (2-5)  the  zeros  and  poles  are  all  complex 
conjugate  for  m  even  and  for  m  odd,  one  of  the  zeros  and  one 
of  the  poles  will  be  negative  real. 

Let  m=3,  the  zero's  and  pole's  pattern  is  plotted  in 
Figure  (2-2) . 
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z-plane 


R  z 
e 


Figure  (2-2)   Pole-zero  pattern. 


where 


z,    is  complex  conjugate  of  z1 , 
T     is  sampling  frequency, 


u 


0 


is  a  notch-frequency  which  the  one  in  which  the  first 
zero  and  pole  occur. 
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For  m=3  it  is  given  by 


'1L 

zl=   ^  e 


•  IT 

Pi  =  S*x      e 


After  multiplying  the  numerator  and  the  denominator  of  (2-3) 
by  z   and  factoring,  the  following  zero-pole  format  is 
obtained 

(z-z,)  (z-z,  )  ...  (z-z?p)  (z-z?p  ) 

H(z)  =  ± i, 2-± 2±^  (2-6a) 

(z-p1)  (z-p1  )  ...  (z-P2£)  (z-P2£  ) 

where  m  is  even  and  given  by 

m  =  21 

I   =    1,  2,  . . . ,  N. 
If  m  is  odd,  called  m  =  2&+1,  the  format  becomes 

(z-z,)  (z-z-.  )  ...  (z-z?p)  (z-z?o  )  (z-z95  -.) 
H(z)  =  ± ^ ^ ^ ZJ6  +  1_   (2_6b) 

(z-Pl)  (z,Pl  )  ...  (z-p2£)  (z-p2£  )  (z-p2£+1) 

where  z?p  +  ,  and  P?p+-i  are  negative  real. 

The  condition  for  the  system  being  stable  is: 


E  |H  (ejU)  T)|  <  -  (2-7) 

o)=0 
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and  the  steady-state  response  of  any  stable  system  with 
transfer  function  (2-6)  to  the  sinusoidal  input, 
u(K)  =  Sin  KooT ,  is  given  by 


where 


y(K)  =  |H  (ejKajT)  |  Sin  (KooT  +  9) 


and 


H  (eja)T)  -  H  (z)|  z  .  eJo)T 


0  =  angle  {H  (e^wT)}  . 


•    rp 

The  term  |H  (eJ   )|  is  called  the  system  gain  factor  since 
it  characterizes  the  sinusoidal  gain  response  property  of 
system  (2-6).   Hence,  the  steady  state  sinusoidal  response 
is  readily  obtained  by  evaluating  the  system  transfer  func- 
tion  at  z  =  e^   ,  where  oj  is  the  radian  frequency  of  the  input 
sinusoid  and  T  is  the  sampling  period.   One  can  determine 
the  system's  gain  and  phase  angle  factor  in  a  straight- 
forward manner  by  using  a  geometrical  approach.   This  approach 
gives  the  advantage  of  relationship  between  the  locations 
of  zeros  and  poles  of  the  system  transfer  function  and  its 
gain-factor  behavior.   In  the  design  of  digital  filter,  it 
will  enable  us  to  design  a  system  whose  gain-factor  approxi- 
mates any  desired  sinusoidal  discrimination  behavior. 
Substituting  z  =  e-'a)   into  (2-6b), 
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If! 

O 


§ 


V? 

— « 


o 

— * 
C) 


§ 


§_ 


OjC 


0-Q 


^T 


02- 


32j 


K-SCRLE=5.00E-01  UNITS  INCH- 
y-SCHLE«S.00E-01  UNITS  INCH- 


Figure  (2-5)   Frequency  magnitude  plot  for  m=3 
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^e    "zi)(e    "zi  }---(e  „;z2^)(e  „"z2^)(e  ,„    &  3 

(eJ   -Pl)(eJ   -Pl  )...(eJ   -p^)  (eJ   -p2£  )  (eJ   "P2£  +  1) 


HCe^1)  = 


where 


(2-8) 


m 


/5-  eJ(2n*l) 

m 


m 


1       1  ^   *^  j    •  •  •  i   *..  *■■ 

n  =  0,  1,  . . . ,  1% 

m  =  21   +  1 


m 


2&  +  1 


m 


Similarly- 


let 


m, —   j(2n+l)- 
p .  =   /a   e J  ^    '  m 


p2£  +  l 


m 


/T" 

m 


a  -  =  e 


jwT  _  m/5-  j(2n+l) 
m 


TT 


0  iwT   my —  i  (2n  +  l) 

$.=eJ    -   /a   eJ  ^    ^m 
i  m 


a2£+l  "  e     +   ^ 
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c        jo)T   m  / — 


Hence  (2-8)  becomes 


H(ejwT) 


a. 


6. 


a 


6. 


a 


2£ 


a 


2% 


3 


2£ 


8   * 

p2£ 


a 


2£  +  l 


(2-9) 


3 


2£  +  l 


It  is  obviously  the  system  gain  factor  at  radian 
frequency  go  is  seen  to  depend  on  the  distance  by  which  the 
system  transfer  function's  zeros  and  poles  are  located  from 
the  critical  point  e        .      Thus,  if  one  wishes  to  completely 
reject  the  sinusoid  sin  IcgoT  ,  it  is  necessary  that  the  trans- 
fer function  have  the  radian  frequency  go  increases  from 
zero  to  infinitive,  one  sees  that  the  zero  of  transfer 
function  (2-3)  will  fall  at  the  first  of  the  zeros  of 
[(2n+l)— ]  radians  for  n=0  and  repeats  for  n=l,2,...  as  in 
Fig.  (2-3). 


Figure  (2-3)   Gain  frequency  behavior  of 
digital  notch-filter. 
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By  considering  pole-zero  pattern  in  Fig.  (2-2)  and  the 
equation  (2-9)  one  wishes  to  reject  completely  at 

z  =  e^     m,  the  distance  a.    must  be  equal  to  zero. 
One  knows  that 


m/rr-     j(2n+l)- 
1  m 


hence 


_   i(2n+l)-   myr—  j(2n+l)- 
a .  =  eJ  ^    -'m  -   /b  eJ  K  Jm 

l  m 


From  the  above  equation, 


a.  =  0      implies  b   =  1 
l  ^      m 


and  for  the  gain  of  (2-9)  is  constant  for  all  frequencies 

except  one  at  z  =  e-*  *■  n   -^m  .   One  should  choose  the  coeffi- 
cient, a  ,  as  close  to  unity  as  possible:  the  closer  a   is 
to  unity  the  closer  the  notch-filter  approaches  the  ideal 
notch-filter  as  shown  in  Fig.  (2-1). 

Looking  back  in  Fig.  (2-2),  the  parameter  A  =  1-a  ,  which 

controls  the  notch-bandwidth,  when  a   close  to  one, 

m 

A  close  to  zero  and  so  does  the  notch-bandwidth.   Ideally, 
it  is  desired  to  make  A  as  small  as  possible,  but  in  prac- 
tice considerations  one  is  not  able  to  make  A  arbitrarily 
small  (A  will  be  discussed  in  section  D  of  this  chapter) . 

The  above  considerations  can  be  made  clear  by  analyzing 
the  transfer  function  (2-3). 
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C.   ANALYSIS  OF  MAGNITUDE  RESPONSE 

It  is  discussed  here  the  magnitude  response  as  a 

function  of  frequency  and  the  coefficients  a   and  b  . 

n    }  mm 

The  equation  (2-3)  is  rewritten  as  follows: 

z   +  b 
H(z)  =  — ^  (2-10) 

z   +  am 

T     _  jwT  and  substitute  into  (2-10),  it  is  given  by 

lc t  z  —  e 

.  _    eJmuT  +  bm 

m 

i  x 
Using  eJ      =  cosx  +  jsinx,  one  gets 

-  ™     (Cos  mwT  +  b  )  +  j  sin  nuoT 
H(e:,a,T)  =  - (2-12) 

(Cos  mwT  +  a  )  +  i  sin  mwT 

v  m    J 

The  squared  magnitude  of  H(e3   )  will  be 

2  2 

.    „      „         (Cos    mooT   +   b    )      +    sin      mooT 

|H(eJ       )|       =   2 2 

(Cos   mcoT   +    a    )      +    sin      mcoT 


or 


2 

-    rp      ~        1    +   b        +    2b    Cos   mwT 

H(e3ajT)|2    =    ^ (2-13) 

1    +   a        +    2a   Cos   mwT 
m 


Take    derivative   of    (2-13)    with   respect    to   w   and    let    the 
numerator   equal    to    zero 


19 


or 


2mb  T  sin  mwT  [1  +  a    +  2a   Cos  mwT] 
m  L     m      m        J 


+  2ma  T  sin  mcoT  [1  +  b  2  +  2b   Cos  mcoT]  =  0 
m  L     m      m        J 


f-2mb  T) (1+a  2)  sin  mu)T  +  2ma  T  (1+b  2)  sin  mwT  =  0 

v    m      m  y  mm' 


Hence 


sin  mwT  [a  (1  +  b  2)  -  b  (1  -  a  2)]  =  0     (2-14) 
L  m      m      m      m  J  J 


Suppose  the  term  in  the  bracket  is  not  equal  to  zero. 
The  frequency  at  which  the  squared  magnitude  is  maximum  or 
minimum  will  satisfy  the  equation  (2-14) 


niT 
sin  mcoT  =  0      implies     coT  =  — 

v  m 


n=  0  , 1 ,  2  ,  .  .  . 


At  n=0,  one  has  ooT=0,  the  magnitude  has  the  maximum 

value  and  at  n=l,  ojT=—  ,  the  magnitude  has  the  minimum  value 

'     m  '       b 

7T 

and  co=—  is  defined  as  a  notch-  frequency ;  when  n=2,3,.... 

the  magnitude  function  just  repeats  as  in  Fig.  (2-4). 

If  one  wishes  to  get  the  constant  gain  except  at  frequency 

77 

w=n— j    >  n-1,3,5, . . . ,  as  in  Fig.  (2-3),  the  coefficients  a 

and  b   are  selected  such  that  this  result  is  achieved, 
m 

To  do  that  in  equation  (2-14),  one  lets  the  factor  in  the 
bracket  always  equal  zero,  i.e., 
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Gain 


Figure  (2-4)   Gain  frequency  behavior  with  some 

values  of  a   and  b  . 
m      m 


a   (1+b  2)  -  b   (1+a  2)  =  0 


or 


2  2 

a   +  a  b    -  b   -  b  a   =  0 

m    mm     m    mm 


a  -b   =  a  b   (a  -b  ) 

mm    m  m  *•  m  nr 


(2-15) 


Hence 


a  b   =  1 

m  m 


We  know  that  the  system  is  stable  when  all  the  roots 
of  characteristic  polynomial  or  transfer  function  of 
system  are  inside  the  unit  circle  of  z-plane. 

By  that  condition,  the  coefficient  a   must  be 

J  '  m 


m 


As  stated  in  section  B,  the  coefficient  b  has  the  value  of 

m 

one  for  completely  rejecting  at  the  notch-frequency.   One 
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sees  that  with  b  =1  and  a   <  1,  the  equality  (2-15)  cannot 
be  satisfied.   But  we  can  approximate  (2-15)  by  letting  the 
coefficient  a  very  close  to  one. 

One  sees  that  both  statements  by  analyzing  the  transfer 
function  (2-3)  and  considering  the  pole-zero  pattern  in 
Fig.  (2-2)  are  appropriate. 

The  transfer  function  (2-3)  now  becomes: 

i  +   "m 
H(z)  =  i 5-^  (2-16) 

1  +  az 

where  the  subscript  m  is  deleted  and  its  squared  magnitude 
is  given  by: 


!„,  jo)T.  i2    2  +  2  Cos  mwT  fn    -,  ^^ 

|H(eJ   )  |   =  ~ (2-loa) 

1  +  a   +  2a  Cos  mwT 

The  plot  of  |H(e-*   )|  for  particular  value  of  a  vs  radian 
frequency  to  for  m=3  in  Fig.  (2-5).   The  half  power  frequen- 
cies are  that  at  which  the  gain  is  3dB  down  or 

|H(e^T)|  -  \ 


or 

2  +  2  Cos  mojT =   1 

1  +  a2  +  2a  Cos  mwT     2 


(2-17) 


Solve  (2-17)  for  go 


2 
4  +  4  Cos  mooT  =  1  +  a   +  2a  Cos  mwT 


(4-2a)  Cos  mooT  =  a   -  3 
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Hence 


Cos  mojT  = 


4  -  2a 


(2-18) 


There  are  two  solutions  go,T  and  w?T  which  are  symmetric 

to  wT=—  ,  the  notch- frequency  of  filter,  as  illustrated  in 
Fig.  (2-5) 


1 

Gain 

1 

707 

r- 



\  B  7 

0J,T     7T_     00?T 

m 


wT 


Figure  (2-6) 

where  B  is  a  notch-bandwidth. 

The  notch-bandwidth  B  will  be 


■x—  =   T(u)-    -    oo,  )    =    2    [—  -    —  Cos         (-. o— )  ] 

2tt  ^2  1  Lm        m  M    -    2a ' J 


or 


d  47T     r 

B     =     I  TT 

m       L 


r     -1    i-a 

Cos         ( 


4    -    2a 


)] 


(2-19) 
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If  one  inserts  the  parameter  A=l-a  into  equation  (2-19), 
the  bandwidth  will  be  a  function  of  A 

B  =  11  [,  -  Cos"1  j1  'J]2    \h 
m   L  4-2(1  -  A) J 


or 


„  _-l  A2    2A 


B  "  m-  C*  "  C0S       2A  ;  2   1         C2-20) 

From  (2-19)  and  (2-20),  the  limit  of  the  bandwidth  B 
when  the  coefficient,  a,  approaches  one  (or  A  approaches 
zero)  is  zero.   This  is  valid  in  theory;  in  fact,  there  is 
some  limit  for  the  coefficient,  a,  (also  for  A)  in  practical 
view  (this  phenomenon  will  be  seen  in  Chapter  V) . 

D.   LOWER  BOUND  OF  A  DUE  TO  FINITE  PRECISION  ARITHMETIC 

Because  of  the  limitations  caused  by  finite  wordlength 
in  all  current  digital  computers,  a  real  number  must  be 
represented  by  using  a  finite  number  of  bits  in  either  the 
fixed-point  or  floating-point  form.   The  error  introduced 
by  such  a  representation  is  discussed  in  this  section.   Only 
binary  arithmetic  will  be  discussed  here. 

In  a  fixed-point  representation,  assume  a  number  x  has 
been  normalized  such  that  |x|  <_   1  can  be  represented  in  twos 
complement  form  as 

oo 

x  =  b.  +   E   b. 2"1 

where 

b.=l  or  0  and  b„  is  the  sign  bit. 
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To  approximate  x  by  a  t-bit  register  length,  one  of  two 
processes  is  generally  utilized:   rounding  or  chopping. 

In  rounding,  a  1  or  a  0  is  first  added  to  the  t-th  bit, 
(x)t_-i>  according  to  whether  the  (t+l)-th  bit,  (x)  ,  is 
1  or  0.   Then  only  the  first  t  bits  of  the  result  are  kept. 
In  chopping,  those  bits  beyond  the  most  significant  t  bits 
are  simply  dropped.   Since  the  error  introduced  by  chopping 
is  larger  than  that  introduced  by  rounding,  chopping  arith- 
metic is  not  commonly  used.   Therefore,  the  discussion 
below  is  limited  to  the  rounding  arithmetic  process. 

Let  (x)   be  the  t-bit  representation  of  the  number  x. 
The  error  which  will  be  introduced  is  defined  by 


E  =   (X)t  -  X 

where  |e|  <  2_t  (2-21a) 


It  is  known  that  if  two  t-bit  fixed-point  numbers  are 
added,  their  sum  still  contains  t  bits,  provided  that  there 
is  not  overflow.   Therefore,  under  this  assumption  of  no 
overflow,  fixed-point  addition  causes  no  error.   On  the 
other  hand,  the  product  of  two  t-bit  numbers  may  have  more 
than  t-bits.   Thus,  rounding  is  necessary  if  only  t  bits 
are  to  be  kept  for  future  processing. 

In  the  case  of  a  floating-point  number,  x  is  represented 
as 

x  =  f2e 
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where  f,  a  fraction  between  y  and  1,  is  called  the  mantissa 
and  e  is  a  binary  integer  called  the  exponent.   Let  x  denote 
the  t-bit  mantissa  floating-point  approximation  of  the 
number  x.   The  difference  between  x  and  fl(x)  =  (x)   is  called 
the  round-off  error  and  depends  on  the  size  of  x.   It  is 
therefore  the  best  measure  relative  to  x.   Then  it  is  clear 
that 

(x)   =  x(l  +  e)  (2-21b) 

where  the  relative  error  e  is  bounded  by  -  2   <^  e  £  2 
Unlike  the  fixed-point  case,  both  addition  and  multiplica- 
tion in  floating-point  can  introduce  round-off  error. 

From  the  preceding  section,  it  is  apparent  that  a 
digital  notch  filter  with  notch  at  u)n=— ~r  can  be  implemented 
by  means  of  a  filter  whose  transfer  function  has  a  zero 
located  at  e^O  .   In  order  that  the  gain  of  such  a  filter 
in  the  frequency  domain  is  near  one  for  frequencies  suffi- 
ciently close  to  oon,  it  is  necessary  that  there  exist 
properly  located  poles  which  tend  to  cancel  the  effects 
of  these  zeros.   This  can  be  accomplished  by  selecting  an 
appropriate  A;  from  equation  (2-20)  it  is  obvious  that  the 
notch  of  the  filter  is  directly  dependent  on  the  parameter 
A.   Ideally,  it  is  desirable  that  A  be  made  as  small  as 
possible  in  order  to  make  the  notch  width  small. 

There  are  practical  considerations,  however,  which  when 
taken  into  account,  will  not  enable  the  filter  designer  to 
make  A  arbitrarily  small  (see  Chapter  V) .   In  the  design 
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of  a  notch-filter  of  (2-16),  one  wishes  the  gain  factor 

at  w=— ~r  to  be  zero.   For  the  critical  point  e^w   sufficiently 


mT 


J 


m 


outside  the  circle  of  center  z,  =  e   and  of  radius  A, 
Fig.  (2-6),  the  gain  factor  for  all  practical  purposes  is 
one  . 


Center 


Figure  (2-6)   A  pair  of  first  pole  and  zero. 

This  follows  since,  for  z  =  e-JW   outside  this  circle,  by 
selecting  A  small  enough,  the  ratio  of  the  product  of  the 
distances  from  zeros  and  poles  to  the  critical  point  eJ 
is  effectively  a  constant.   Since  all  digital  computers 
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are  limited  by  finite  wo-rdlengths ,  the  desired  coefficients 
will  be  truncated  and  the  coefficient  inaccuracies  will 
result  as  stated  at  the  beginning  of  the  D  of  Chapter  II. 

The  characteristic  difference  equation  of  a  digital 
filter  (2-16)  is  given  by 

v(K)  =  u(K)  +  u(K-m)  -  av(K-m)  (2-22) 

When  one  wishes  to  implement  (2-22)  on  a  digital  computer, 
the  coefficient,  a,  will  be  stored  with  error  due  to  trunca- 
tion as  (a+e) ,  where  e    is  the  round-off  error  in  (2-21a,b). 
This  will  rise  to  a  shift  in  the  actual  pole  locations  of 
the  transfer  function. 

For  system  stability,  the  condition  placed  on  (a+e)  is: 

a  +  e  <  1 
or  1-a  >  £ 

hence         A  >  e  (2-23a) 

If  one  uses  a  t-bit  wordlength,  (2-23)  will  become: 

A  >  2_t  (2-23b) 

In  the  implementation  of  digital  notch-filter  (2-22) 
with  t-bit  wordlength,  one  must  select  the  parameter 
A  >  2    to  guarantee  that  the  notch-width  is  sufficiently 
narrow  and  to  avoid  the  possibility  of  unbounded  system 
growth.   It  is  noted  that  if  one  counts  on  guard  bit,  the 
wordlength  will  be  (t+1)  bits.   Given  notch-uidth  B,  one 
can  calculate  A  from  (2-20)  and  therefore  t. 
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III.   STATE  SPACE  REPRESENTATIONS  AND  REALIZATIONS 

The  transfer  function  (2-16)  of  an  m-th  order  digital 
filter  is  the  relationship  between  one  input  and  one  output. 
Only  state  and  output  equations  with  single  input  and  single 
output  will  be  considered.   In  the  time  invariant  discrete 
system,  the  transfer  function  (2-16)  can  be  described  by 
a  state  variable  representation  of  the  form 

x(K+l)  =  A  x(K)  +  B  u(K)  (3-la) 

V(K)  =  C  x(K)  +  6  u(K)  (3-lb) 

where 

A  is  an  m  x  m  matrix, 

B  is  an  m  x  1  matrix 

C  is  a  1  x  m  matrix 

6  is  a  1  x  1  matrix  or  scalar. 
Taking  the  z-transform  of  (3-la)  and  solving  for  X(z)  yields 

X(z)  =  [I  -  Az"1]"1  B  U(z)  (3-2) 

where  I  is  an  m-dimension  identity  matrix.   Then,  taking 
the  z-transform  of  equation  (3-lb)  and  substituting  (3-2) 
for  X(z)  yields 

V(z)  =  (C  [I  -  Az"1]"1  B  +  6}  U(z) 

=  H(z)  U(z)  (3-3) 
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where 


H(z)  =  C  [I  -  Az"1]"1  B  +  6 


(3-4) 


H(z)  is  the  transfer  function  of  system  (3-1).   The 
discrete  system  is  considered,  as  stated  above,  for  only 
one  input  and  one  output.   Hence,  the  matrix  parameters 
A,  B,  C  and  5  have  the  forms: 


m  x  m 
A  = 


m  x  1 

B  = 


11 


ml 


lm 


m  x  m 


(3-5) 


1  x  m 
C  = 


[C. 


•  •  C  J 

m 


6  =  d 
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There  are  many  ways  to  realize  the  digital  filter  (3-1) 
Some  of  them  are  the  canonic  realizations.   The  term 
"canonic"  is  used  in  the  sense  of  minimizing  the  number  of 
operations  required.   This  minimization  implies  a  minimum 
number  of  coefficients  of  the  difference  equations,  which 
is  the  source  of  noise. 

If  one  wants  to  describe  the  system  (3-1)  by  a  matrix 
form  as : 


S  = 


B 


C   i   6 


(3-7) 


where  S  is  called  characteristic  matrix  system. 

The  matrix  parameters  A,  B,  C  and  5  are  respectively 
given  as  follows: 


A  = 


0  - 


(m-1)  x  (m-1) 


a 
0 


B  = 


(3-6) 
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c  =  [ill  o  .  .  .  0] 


«-i 


To  describe  a  system  in  the  matrix  form,  its  state 
equations  will  have  the  form: 


x1  (K+l)  =  -  a  xm  (K)  +  u  (K) 


x2  (K+l)  =    x±    (K) 


(3-7) 


xm  (K+l)  =  x   -  (K) 

m     J  m-1  v  J 


v  (K)     =  (5ll)  X;L  (K)  +  \   u  (K) 


The  corresponding  realization  of  (3-7)  has  the  form 
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1 
a 

i(K) 

, 

,} 

c 

rr\ 

a-1 
a 

J. 

•  ^    v,c 

X 

_,         r  —         -^ 

V 

V 

Figure  (3-1)   Non-canonic  realization  of  m-th  order 

digital  notch-filter. 

The  realization  in  Fig.  (3-1)  is  not  canonic;  the  purpose 
of  this  paper  is  to  use  only  one  coefficient.   There  are  many 
disadvantages  in  the  above  realization.   Among  these  are 
the  production  of  more  noise,  higher  cost  and  greater  network 
complexity.   Hence,  the  matrix  form  (3-7)  is  a  poor  way  to 
describe  the  system.   This  leads  to  the  consideration  of 
a  search  technique  for  the  canonic  realization  as  stated 
in  section  A.   From  this  concept,  one  can  realize  the  digital 
notch-filter  (2-6)  as  follows: 


33 


u(K) 


v(K) 


Figure  (3-2)   Canonic  realization  of  digital 

notch-filter. 

and  the  corresponding  state  equations  of  realization  (3-1) 
are  listed  below 


xx  (K)     =  -a  xm  (K+l)  +  u(K) 
Xj  (K+l)   =   xx  (K) 


x2  (K+l)   =  x±    (K+l) 


(3-8) 


xm  (K+l)   =  xm  ,  (K+l) 

m     ' 

v  (K) 


=  xm  (K+l)  +  Xl(K) 


One  sees  the  realization  in  Fig.  (3-2);  with  this 
technique  there  is  only  one  coefficient  and  it  provides 
the  minimum  number  of  the  coefficients  of  the  m-th  order 
recursive  digital  filter  one  can  achieve.   Therefore,  in 
noise  analysis,  if  one  consider  only  the  source  of  noise 
produced  by  the  multiplications  of  the  coefficients,  it 
is  not  necessary  to  assume  the  sources  of  noise  are 
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uncorrelated,  since  there  is  only  one  source  of  noise  in 
the  realization.   By  this  advantage,  in  the  design  of  high 
order  of  digital  notch-filter  it  is  not  necessary  to  imple- 
ment by  cascading  first  or  second  order  sections  as  is 
often  done  at  present. 

The  "dual"  canonic  realization  of  Fig.  (3-2)  can  be 
derived.   The  so-called  "transposed  realization"  has  been 
shown  by  Jackson  [Ref.  4]  to  be  given  in  the  form 


*v(K) 


Figure  (3-3)   Transposed  canonic  realization 
of  digital  notch-filter. 

The  corresponding  equations  are  given  by: 
xx  (K)    =  -a  u(K)  +  u(K) 
x±    (K+l)  =   X;L  (K) 
x2  (K+l)  =   X;L  (K+l) 


(3-9) 


x   (K+l)  =   x   ,  (K+l) 

m    K        '  m-1  v    J 

v(K)      =   xm  (K+l)  +  u(K) 
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In  the  state  equations  (3-8)  and  (3-9),  the  value  of 
the  coefficient,  a,  is  taken  between  0  and  1.   It  is  known 
that  the  notch-gain  and  notch-width  of  a  digital  notch- 
filter  are  functions  of  the  coefficient  a.   The  question 
which  arises  here  is:   what  is  the  optimum  value  of  the 
coefficient  a  which  satisfies  the  desired  characteristics 
of  the  notch-filters  and  how  large  is  the  ripple  in  their 
passband? 

These  questions  are  investigated  in  the  following 
chapters  by  simulating  the  digital  notch-filter  with  the 
digital  computer  IBM  360/70.   One  can  see  the  difference 
between  theory  and  practice  in  Chapter  III. 

For  example,  the  gain  in  equation  (2-9)  at  notch 
frequency,  in  theory,  is  a  decreasing  function  when  the 
coefficient  a  is  increasing.   However,  in  practice,  some 
values  of  coefficient  a  make  the  notch-gain  change  non- 
monotonically  because  of  round-off  noise  in  iteration 
computations  in  the  program. 
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IV.   IMPLEMENTATION 

A.   COMPUTER  IMPLEMENTATION 

In  this  chapter,  a  digital  notch-filter  which  is 
characterized  by  the  state  equations  (3-8)  is  implemented 
on  the  IBM  360/70.   Plots  of  the  gain-frequency  response 
and  the  time  response  are  provided. 

The  input  signal  is  a  sine-wave,  i.e.,  a  signal  of 
single  frequency,  which  sampled  every  T  seconds.   The 
corresponding  gain  for  each  input  sine-wave  is  computed 
when  the  steady  state  is  reached,  by  the  ratio  of  the  value 
of  the  peak  input  sine-wave  divided  by  the  value  of  peak 
output  response.   Then  the  input  sine-wave  is  scanned  over 
the  frequency  range  of  the  notch-filter  and  the  gain- 
frequency  response  is  plotted.   It  is  noted  that  the  notch 
of  the  filter  falls  at  the  first  zero  of  filter;  that  is, 

-TT 

at  z-,  =  e   .   If  one  desires  the  notch-filter  at  oon 
radians,  the  sampling  time  T  is  computed  as  follows 

^0T    Jm 
e     =  e 

or 

0    m 
Hence 

TT 


T  = 


rnw^. 
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or 


T  = 


2mi 


0 


(4-1) 


to 


0  . 


where  fQ  =  ■=—  is  the  notch- frequency  in  hertz. 

The  following  sections  contain  implementations  of 
second,  third,  and  sixth  order  digital  notch-filters  on 
the  IBM  360  computer  using  single  precision  arithmetic. 

1 .   Second  Order  Digital  Notch-Filter 

The  results  of  m-th  order  digital  notch-filter 
in  Chapters  II  and  III  are  applied.   The  second  order 
digital  notch-filter  is  characterized  by  the  transfer 
function 


H(z)  = 


1  +  z 


-2 


1  +  az 


(4-2) 


and  it  is  realized  as  follows 


^(K+l) 


x2(K+l) 


v(K) 


Figure  (4-1)   Canonic  realization  of  second  order 
digital  notch-filter. 
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The  corresponding  state  equations  are  given  by 

x-l  (K)    =  -a  x2  (K+l)  +  u  (K) 
x1    (K+l)  =  x±    (K) 


x2  (K+l)  =   xx  (K+l) 

V  (K)     =   x2  (K+l)  +  Xl  (K) 


(4-3) 


If  one  wishes  a  notch  frequency  at  fn  =  60  hz,  the 
equation  (4-1)  gives  the  sampling  frequency  T  as  follows 

T  =  -^r-  =  9  -  I   -  ,n   =    4.166666667  x  10"3 
zmf^    2  x  2  x  60 

The  figures  from  (4-2)  to  (4-10)  illustrate  the 
gain  frequency  responses  in  dB  with  different  coefficients 
of  one  thousand  iterations.   It  is  noted  that  at  one  thousand 
iterations  the  steady  state  of  frequency  response  is  not 
yet  reached.   These  plots  merely  depict  the  behavior  of 
the  gain-frequency  response  with  changing  coefficients. 

The  result  of  the  simulating  program  gives  the  steady 
state  of  frequency  response  at  about  six  thousand  iterations. 
They  are  shown  in  Figs.  (4-2)  to  (4-17)  corresponding  to 
some  typical  coefficients.   And  the  notch-gains,  notch- 
widths,  and  maximum  passband-ripple  in  dB  of  the  steady 
state  gain- frequency  response  are  tabulated  in  Table  (4-1). 
Fig.  (4-11)  shows  the  plot  of  notch-width  vs.  coefficient 
in  theory  by  evaluating  the  equation  (2-19). 
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The  time  response  of  a  second  order  digital  notch- 
filter  when  the  input  is  a  sinusoidal  signal  of  frequency 
60  Hz,  i.e.,  at  the  notch  frequency  given  below.   In  theory, 
the  notch-filter  with  60  Hz  notch  will  immediately  reject 
the  above  input  signal,  but  in  practice  the  time  response 
sinusoidal  signal  will  be  damped  in  some  finite  time,  small 
or  large  according  to  the  coefficient  and  order  of  notch- 
filters  one  utilizes.   Figures  (4-18)  to  (4-21)  show  the 
time  response  of  second  order  60  Hz  notch-filters  with  60  Hz 
sinusoidal  input  signal  applied. 
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Figure  (4-2)   Frequency  response  of  second  order 
digital  notch-filter.  a=0. 90714359 
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Figure  (4-4)   Second  order  digital  notch-filter 
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Figure  (4-5)   Second  order  digital  notch-filter 
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Figure  (4-6)   Second  order 

digital  notch-filter. 
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Figure  (4-7)   Second  order  digital  notch-filter. 

a=0.9815 


46 


§ 


2Z2 


Q2: 


r^r 


!J 


0..1 


^iCl 


9 


9 


84.27  dB 


K-SCflLE^2.QQE+Ql.  UNITS  INCH. 
^-SCRL£>-2.00rc+Ql  UNITS  INCH. 


Hz 


Figure  (4-8)   Second  order  digital  notch-filter. 
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Figure  (4-9)   Second  order  digital  notch-filter 
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Figuer  (4-10)   Second  order  digital  notch-filter 
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Figure  (4-11)   Plot  of  notch-width  vs  coefficient, 

in  theory, of  equation  (2-19). 
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Figure  (4-13)   Steady-state  gain  frequency  response 
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Figure  (4-14)   Steady-state  frequency  response 

of  second  order  notch-filter. 
a=0.99 
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Figure  (4-15)   Steady-state  gain  frequency  response 

of  second  order  notch-filter. 
a=0.995 
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Figure  (4-16)    Steady-state  gain  frequency  response 

of  second  order  notch- filter . 
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Figure  (4-17)   Steady-state  gain  frequency  response 
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It  is  noted  that  the  passband  ripple  is  large  for 
frequencies  very  close  to  the  notch- frequency .   Further, 
about  one  hertz  from  the  notch,  all  passband  starts  less 
than  0.09  dB  in  magnitude  and  smoother  for  going  farther 
the  notch. 

2.   Third  Order  Digital  Notch-Filter 

The  third  order  digital  notch-filter  is  charac- 
terized by  the  transfer  function 


H(z)  = 


1  +  z 


1  +  az 


-3 


(4-4) 


and  realized  as  follows 


u(K) 


x-^K) 


x1(K+l)  rv   x2(K+l)  jv   x3(K+l) 


-a 

i i 


v(K) 


Figure  (4-22)   Canonic  realization  of  third-order 
digital  notch-filter. 
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The  corresponding  state  equations  are  derived  as 

xx  (K)    =  -  a  x2  (K+l)  +  u  (K) 

xx    (K+l)  =  x±    (K) 

x2  (K+l)  =  xx  (K+l)                  (4-5) 

x3  (K+l)  =  x2  (K+l) 

v  (K)     =  x3  (K+l)  +  X]L  (K) 

and  the  sampling  time  required  for  60  Hz-notch  is  given  as 


=  2.77777778  x  10"3  sec. 


2mfQ    2  x  3  x  60 


Similar  to  the  results  of  the  second  order  case, 
the  frequency  responses  of  a  third  order  digital  notch- 
filter  are  shown  in  Figs.  (4-23)  to  (4-31)  with  one  thousand 
iterations.   The  steady  state  of  third  order  is  reached  in 
about  ten  thousand  iterations.   The  steady-state  gain 
frequency  responses  are  illustrated  in  Figs.  (4-32)  through 
(4-37). 

The  time  responses  of  third  order  digital  notch- 
filter  are  shown  in  Figs.  (4-38)  to  (4-41).   Similarly, 
Table  (4-1)  of  second  order,  Table  (4-2)  of  third  order 
is  tabulated. 
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Figure  (4-24)   Frequency  response  of  third  order 

notch-filtre . 
a=0. 96901345 
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Figure  (4-25)   Frequency  response  of  third  order 

notch-filter. 
a=0.979 
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Figure  (4-26)   Frequency  response  of  third  order 

notch-filter. 
a=0.98 
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Figure  (4-27)   Frequency  response  of  third  order 

notch-filter. 
a=0.981 
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Figure  (4-28)   Frequency  response  of  third  order 

notch-filter. 
a=0.9815 


^9 


n 

O 


-i 

(  > 

I 


O 
I 


ft 
r> 

Y 


o 
o 

1 


QU 


~*\fr 


Hz 


r> 


47.76    dB 


k-scrle: --.-2.00E+01  units  inch, 

V-SCRLf. --.-I.  D0E+D1  UNITS  INCH. 


Figure  (4-29)   Frequency  response  of  third  order 

notch-filter. 
a=0.982 
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Figure  (4-30)   Frequency  response  of  third  order 

notch-filter. 
a=0. 98298407 
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Figure  (4-31)   Frequency  response  of  third  order 

notch- filter. 
a=0.  99845153 
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Figure  (4-33)   Steady-state  frequency  response  of 

third  order  notch-filter.  a=0.985 
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Figure  (4-35) 


Steady-state  frequency  response  of 
third  order  notch- filter . 
a=0.995 
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Figure  (4-36) 


Staedy-state  frequency  response  of 
third  order  notch-filter. 
a=0.999 
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Figure  (4-37)   Staedy-state  frequency  response 

of  third  order  notch-filter. 
a=0.9999 
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3.   Sixth  Order  Digital  Notch-Filter 

The  sixth  order  digital  notch-filter  is 
characterized  by  the  transfer  function 


1  +  z-6 
H(z)  =  i Z-^  (4-6) 

1  +  az 


One  often  realizes  the  sixth  order  digital  filter 
by  cascading  three  sections  of  second  order  filters.   This 
is  done  to  reduce  the  sensitivity  of  pole  positions  due 
to  truncation  because  of  round-off  error  of  iteration 
computations  in  the  simulation  program. 

The  sixth  order  digital  notch-filter  characterized 
by  transfer  function  (4-6)  is  discussed  here.   There  is 
only  one  source  of  error;   the  noise  produced  in  the 
computer  implementation  is  small  compared  with  the  general 
form  of  sixth  order  digital  notch-filter  of  thirteen 
coefficients.   Therefore,  one  can  implement  the  sixth 
order  digital  notch-filter  (4-6)  by  direct  canonic  reali- 
zation as  given  by  Fig.  (4-42). 
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The  corresponding  state  equations  are  given  by: 


a  x,  (K+l)  +  u  (K) 


xl 

(K) 

= 

xl 

(K+l) 

= 

x2 

(K+l) 

= 

x3 

(K+l) 

= 

X4 

(K+l) 

= 

X5 

(K+l) 

= 

X6 

(K+l) 

= 

v(K 

= 

(K) 

(K+l) 

(K+l) 

(K+l) 

(K+l) 

(K+l) 

(K+l)    + 

Xl 

(K) 

(4-7) 


One  still  requires  the  notch-frequency  at  fn  =  60  Hz 
the  equation  (4-1)  gives  the  sampling  frequency  T  as  follows 


T  = 


2mf 


0 


2  x  6  x  60 


=  1.388888889  x  10"3  sec. 


It  is  noted  that  the  sampling  time  for  sixth  order  is  one 
half  of  that  for  the  third  order  case. 

Figures  (4-43)  to  (4-46)  illustrate  the  time  response 
of  a  sixth  order  digital  notch-filter  with  input  sinewave  of 
60  Hz  frequency. 

Figures  (4-47)  to  (4-49)  are  frequency  response  plots 
of  one  thousand,  two  thousand  and  two  thousand  and  five 
hundred  iterations,  respectively,  corresponding  coefficient 
a  =  0.98298407. 

It  is  noted  that  the  sixth  order  digital  notch-filter 
requires  about  fifteen  thousand  iterations  to  reach  the 
steady-state.   The  steady-state  gain  frequency  responses  are 
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illustrated  in  Figs.  (4-50)  to  (4-55)  and  some  numerical 
values  are  tabulated  in  Table  (4-4)  . 
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Figure  (4-47)   Sixth  order  digital  notch-filter 

at  60  Hz  notch  and  1000  iterations 
a=0. 98298407 
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Figure  (4-48)   Sixth  order  digital  notch- filter 

at  60  Hz  notch  and  2000  iterations 
a=0. 98298407 
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Figure  (4-49)   Sixth  order  digital  notch-filter 

at  60  Hz  notch  and  2500  iterations 
a=0.  98298407 
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Figure  (4-50)   Steady-state 
frequency  response  of 
sixth  order  notch-filter 
a=0.95 
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Figure  (4-51)   Steady-state 
frequency  response  of  sixth 
order  notch-filter. 
a=0.985 
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Figure  (4-52) 
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Staedy-state  frequency  response 
of  sixth  order  notch-filter. 
a=0.99 
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Figure  (4-53)   Steady-state  frequency  response 

of  sixth  order  notch-filter. 
a=0.995 
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Figure  (4-54)   Staedy-s tate  frequency  response 

of  sixth  order  notch-filter. 
a=0.999 


^ 


! 


:^a 


3c 


Hz 


2.2    dB 


K-Sr.RLF^2.0CF+0I.  UNITS  INCH.. 
V'-Sr.RLF--5.  OCF.-  DI  UNITS  INCH- 


Figure  (4-55) 


Steady-state  frequency  response 
of  sixth  order  notch-filter. 
a=0.9999 
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One  is  lead  to  attempt  to  realize  the  sixth  order 
digital  notch-filter  by  cascading  three  sections  of  second 
order  as  follows: 


u(K) 

-W-f^Ir 


Figure  (4-56)   Realization  of  sixth  order  digital 
notch-filter  by  cascade. 

The  steady-state  gain  frequency  response  is  plotted 
in  Fig.  (4-57)  and  (4-58)  for  coefficient  a=0.99  and  a=0.995, 
respectively. 

Compare  with  the  results  of  direct  realization  of 
sixth  order  notch-filter  in  Figs.  (4-52)  and  (4-53);  the 
notch-gain  in  realization  by  cascading  is  deeper  than  in 
direct  realization.  But  in  passband  it  does  not  improve 
anything.  It  seems  worse  than  direct  form  by  introducing 
more  noise;  one  can  see  gain  frequency  response  behavior 
for  coefficient  a=0. 9999999  between  them  in  Figs.  (4-59) 
and  (4-60) . 

It  proves  that  the  result  is  contrary  with  the 
statement  of  many  current  lectures:   the  pole-position  of 
direct  form  is  sensible  than  cascade  form  in  realizing 
high  order  digital  filter  [Ref.  1] .   The  result  in 
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Figs.  (4-59)  and  (4-60)  states  that:   the  pole-position 
in  cascade  form  is  much  more  sensible  than  direct  form  in 
realizing  sixth  order  digital  notch-filter  characterized 
by  the  transfer  function  (2-16)  . 
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Figure  (4-57)   Frequency  response  of  sixth  order 

notch-filter  by  cascade. 
a=0.99 
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Figure  (4-58)   Frequency  response  of  sixth  order 

notch-filter  by  cascade. 
a=0.995 
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Figure  (4-59)   Frequency  response  of  sixth  order 

notch-filter  by  direct  form. 
a=0. 9999999 
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Figure  (4-60)   Frequency  response  of  sixth  order 

notch-filter  by  cascade  form. 
a=0. 9999999 
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4 .   Observation  and  Discussion 

In  theory,  by  inspection,  the  magnitude  response 
in  equation  (2-9)  and  the  notch-width  in  equation  (2-19) 
vs.  coefficient  of  m-th  order  digital  notch-filter  charac- 
terized by  transfer  function  (2-16),  one  sees  that  when 
coefficient,  a,  approaches  one,  the  magnitude  of  notch-gain 
and  notch-width  will  be  smaller  [see  Fig.  (4-11)].   From 
the  magnitude  response  plots,  it  is  seen  that  theory  and 
practice  are  in  agreement  whenever  the  steady-state  output 
response  is  reached.   This  can  be  seen  in  the  steady-state 
frequency  responses  in  Figs.  (4-12)  to  (4-17)  of  second 
order,  Figs.  (4-32)  to  (4-37)  of  third  order  and  Figs.  (4-50) 
to  (4-55)  of  sixth  order.   Some  of  numerical  results  of 
steady-state  response  are  tabulated  in  Tables  (4-1),  (4-2) 
and  (4-3)  . 

For  the  case  of  intermediate  state  (e.g.,  state  of 
one  thousand  iterations- -the  steady-state  responses  are  not 
yet  reached),  the  theory  and  practice  appear  to  disagree. 
One  can  see  this  in  frequency  response  plots  in  Figs.  (4-2) 
to  (4-10)  of  second  order.   These  phenomena  are  a  result 
of  the  random  errors  produced  by  truncating  coefficients 
and  round-off  error  in  iteration  computations  in  the  simu- 
lating program  due  to  finite  register  length  of  the  computer. 
In  this  case,  it  is  due  to  limitations  of  single  precision 
arithmetic  on  IBM  360.   A  second  explanation  of  the  apparent 
deviation  from  the  theoretical  stems  from  error  in  computa- 
tion of  notch-gain  in  computer  program: 


108 


r-„  •       on  t      r  I  Peak  of  input  is 

Gam.r,  =  20  log..  n  f  „ — , i £ — -  ) 

AB        &10   'Peak  of  output ]J 


The  peak  of  input  is  always  one,  but  in  the  computer  program, 
the  peak  of  input  and  output  are  picked  simultaneously  and 
it  is  difficult  to  choose  the  sampling  time  such  that  it 
exactly  corresponds  to  input  peak;  this  is  demonstrated  in 
Fig.  (4-56) 


input 
output 


Figure  (4-56) . 

where  the  sampling  time  T  is  fixed  for  prescribed  magnitude 
response  of  desired  notch-filter. 
Suppose : 

the  peak  of  input  ^0.9 

the  peak  of  output  ~_  0.0002 
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The  corresponding  notch-gain  will  be: 

GAIN  =  -  73.06  dB. 

Now  let  the  peak  close  to  one;  suppose: 
the  peak  of  input  ~_  0.99 
the  peak  of  output   0.0001 

GAIN  =  -  79.91  dB. 

The  error  thus  realized  approximately  6  dB  for 
the  notch-gain.   By  the  difficulty  of  picking  exactly  the 
peak  of  input  as  well  as  the  peak  of  output,  the  actual 
notch-gains  are  certainly  somewhat  deeper  than  the  values 
listed  in  Tables  (4-1),  (4-2),  and  (4-3).   Also,  the  actual 
passband  ripples  should  be  smoother. 

It  is  known  that  the  second  order  needs  around 
six  thousand  iterations,  the  third  order  needs  around  ten 
thousand  iterations  and  the  sixth  order  needs  around  fifteen 
thousand  iterations  to  the  steady-state.   It  is  seen  by 
considering  the  time-respond  plots  of  second,  third  and 
sixth  order.   It  is  obviously  true  since  higher  order  in 
direct  realization  requires  larger  number  of  delays. 

It  is  noted  that  the  passband-ripple  is  large 
around  notch- frequency  and  smoother  when  going  further. 

The  results  show  that  the  higher  order  of  digital 
notch-filter  improves  the  passband-ripple  (i.e.,  reduced 
ripple)  and  notch-gain  [i.e.,  deeper  notch.   See  Tables 
(4-2),  (4-2),  and  (4-3)]  and  steady-state  frequency  response 
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plots  in  Figs.  (4-50)  to  (4-55)  of  second,  third  and 
sixth  order,  respectively. 

The  sixth  order  notch-filter  requires  about  four 
minutes  to  plot  one  steady-state  gain  frequency  response 
using  fifteen  thousand  iterations. 

One  will  see  that  the  results  in  this  chapter  will 
give  the  easy  way  to  work  in  optimum  problems  due  to  finite 
register  length  on  computer  in  the  next  chapter. 

B.   CCD  IMPLEMENTATION* 

The  charge-coupled  device  (CCD)  was  invented  at  the 
Bell  Telephone  Laboratory  by  Drs .  Boyle  and  Smith  in  1970. 
This  new  development  has  been  receiving  much  attention 
by  industry  because  of  its  important  properties- -low  power 
consumption,  low  noise,  and  delay  capability.   Many  success- 
ful applications  have  been  investigated  such  as  digital 
signal  processing  especially  in  design  of  recursive  digital 
filter.   More  details  of  CCD  can  be  seen  in  many  current 
lectures.   The  CCD  hardware  implementation  of  recursive 
digital  notch-filter  is  beyond  this  research.   The  following 
is  the  sample  of  block  diagram  of  second  order  CCD  digital 
notch-filter: 


*"Charge-coupled  Device  for  Analog  Signal  Processing 
Recursive  Filters  Study,"  Engineer's  Thesis  by  Varachai 
Imsa-ad,  Naval  Postgraduate  School,  December  1974. 


Ill 


u(K) 


v(K) 


Figure  (4-61)   Block  diagram  of  a  second  order 
CCD  digital  notch- filter . 
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V.   OPTIMUM  COEFFICIENT 

In  this  chapter,  a  practical  method  is  described  for 
choosing  the  coefficients  of  a  digital  filter  to  meet 
arbitrary  specifications  of  the  magnitude  characteristic. 

As  shown  above,  all  special  purpose  computers  have 
finite  wordlengths  which  is  a  major  source  of  noise  in  the 
implementation  of  digital  filters.   This  is  because  of  the 
necessity  to  truncation  coefficient  length  and  also  because 
of  rounding  in  the  arithmetic  operations,  especially  in 
the  iterative  multiplications.   This  kind  of  noise  deter- 
mines the  signal-to-noise  ratio  of  the  system.   In  some 
cases,  the  limitation  of  wordlength  of  coefficients  may 
lead  to  a  system  which  does  not  satisfy  the  original 
conditions . 

Especially,  in  the  design  of  recursive  digital  notch- 
filters,  one  has  to  select  the  coefficients  such  that  they 
satisfy  the  prescribed  characteristic  of  notch-gain,  notch- 
width,  and  the  smallest  passband  ripple. 

The  m-th  order  digital  notch-filter  is  characterized 
by  the  transfer  function  (2-16).   It  is  seen  that  the  notch- 
gain,  the  notch-width,  and  the  passband  ripple  are  functions 
of  a  single  parameter- -the  coefficient  "a"  between  0  and  1. 
It  is  possible  to  construct  a  performance  function  of 
single  parameter  in  the  optimization  problem. 
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A.   SEARCH  PROCEDURE 

In  the  design  of  an  m-th  order  digital  notch-filter, 
if  one  wishes  the  prescribed  magnitude  having  the  follow- 
ing characteristics: 

(i)  Notch-gain  is  at  least  -  60  dB. 
(ii)  Notch-width  is  at  least  one  hertz, 
any  coefficient,  a,  in  )0,1(  satisfies  conditions  (i)  and 
(ii) ,  and  that  value  which  produces  the  smoothest  passband 
is  the  optimum  value. 

In  the  method  of  searching,  first,  one  selects  the 
subinterval  of  )0,1(  such  that  any  coefficient  belonging 
to  this  subinterval  will  satisfy  the  conditions  (i)  and 
(ii) ;  second,  searching  is  accomplished  for  the  minimum 
error  function  of  single  variable  for  one  of  the 
subintervals . 

The  subinterval  can  be  found  by  computer  experiment. 
Suppose  for  sixth  order  notch-filter  one  can  consider 
Table  (4-3);  call  (a,,  a2)  this  subinterval: 


al        a2 
0y//////  //ss/  //////  //////////'I \tfM'///  /A. 


As  stated  above,  any  coefficient,  a,  belonging  to 
interval  (a, ,  a?)  will  satisfy  the  conditions  (i)  and  (ii) 

For  sixth  order,  the  values  of  a,  and  a~  are  given  as 
[see  Table  (4-3)] : 


a1   =  0.98    and    a2  =  0.995  . 
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For  third  order  [see  Table  (4-2]  : 

a±   =  0.98    and    a      =  0.99  . 

It  is  noted  that  the  ideal  (errorless)  magnitude 
frequency  response  is  equal  to  zero  dB  over  the  passband 
of  the  notch-filter.   Due  to  the  limited  wordlength  of 
computers ,  the  actual  magnitude  frequency  response  will 
deviate  from  zero  dB  over  the  passband.   If  one  defines: 

c  (a,  ok)  =  1  -  H  (a,  u/)  (5-1) 

as  the  error  function  corresponding  to  coefficient,  a,  and 
radian  frequency  co-  of  sinusoidal  input  signal.   The  error 
function  e  (a,  w.)  is  defined  in  discrete  domain  which  is 
considered  as  approximates  the  mean-squared  error  in  con- 
tinuous domain  of  interval  between  0  and  oo .  . 


03- 

(a,  u.0  ~£-  f        K  (a>  wi)  "  H  (a>  w-,)|2  dwi    C5-2) 
1    wi   0 


where  Hto  (a,  co.)  is  ideal  response  and  has  magnitude  equal 
to  one . 

Rewrite  (5-2)  as 


0). 

(a,  to.)  =  —  /     1 1  -  H  (a,  u.)\*   dec.       (5-3) 


the  system  is  stable,  e  (a,  to.)  has  to  be  convergent. 
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Let  a),  scan  over  the  frequency  range  of  the  passband  of 
notch-filter,  the  error  performance  or  the  total  error 
function,  is  defined  as 


N    2 
E  (a)  =   I   e   (a ,  oj. ) 

i  =  l         1 

or 

N  2 

E(a)  =   Z   [1  -  H  (a,  to . )  ] L  (5-4) 

i=l  x 


where  a,  £  a  £  a„  and  N  is  a  finite  integer  number. 

The  optimization  problem  turns  out  to  be  the  least 
square  method  of  single  parameter.   The  optimum  coefficient, 
a,  is  one  that  makes  the  error  performance  E(a)  minimum. 
The  following  algorithms  are  described  to  search  for  the 
minimum  of  E(a).   First,  the  interval  (a,,  a  )  is  subdivided 
into  N  subintervals  such  that 


(K-l)(a2-a1)        K(a2-ai) 


aK  e  [al  +  N ■  al  +  N ] 


where  K  =  1,  2,  ...,  N. 


Then  compare  the  values  of  E(a„)  and  ECa^,,.  corresponding 


to  successive  value  of  a„  and  a^+1 ,  where 


K(a  -a  )         (K+l)(a?-a  ) 

a    e  [a,  +  k — —  ,  an  +  ~ — — ] 

K+l      1       N        1         N      J 


If    E(aK)  >  E(aK+1) 
the  next  step  will  compare  E  (a„+^  and  E  (aj,  +  J  and  store 
coefficient  a„+,. 
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Or  if  E  (a„)  <_   E  (aK+,)  ,  the  next  step  will  compare 

E  (a^)  and  E  (aK+?-,  and  store  coefficient  a^- ;  the  execution 
continues  until  the  last  subinterval  of  (a-,,  a2)  is  reached. 

The  final  result  will  be  printed  out  as  the  optimum 
coefficient  corresponding  to  the  smallest  value  of  error 
performance  E  (a)  for  a,  <_   a  <_  a?. 

The  greater  the  value  of  N,  the  higher  accuracy  one 
can  get  but  the  price  is  paid  by  spending  more  computation 
time . 

The  plots  of  error  performance  are  illustrated  in  Fig. 
(5-1)  of  second,  third,  and  sixth  order  digital  notch- 
filters.   One  sees  that  the  error  performance  is  monotoni- 
cally  decreasing  with  increasing  coefficient,  i.e.,  the 
greater  the  coefficient,  the  smoother  the  response  in 
passband.   This  result  can  be  extended  for  many  order  of 
digital  notch-filter  characterized  by  the  transfer  function 
(2-16).   The  optimum  coefficient  now  is  so  easy  to  pick  and 
it  is  the  coefficient  a~. 

for  sixth  order    a.      =  0.995 
for  third  order    a~  =  0.99. 

B.   DISCUSSION  AND  CONCLUSION 

One  sees  that  the  error  performance  of  single  parameter 
is  unimodal.   The  optimum  problem  of  this  kind  of  function 
has  been  accomplished  by  many  methods  well-known  such  as 
Fibonacci  Search,  Golden  Section,  etc.   The  Fibonacci  Search 
is  known  as  the  most  efficient  in  optimization  technique  of 
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o 


sixth  order 
Third  order 

Second  order 


K-SCRLF.:"-2.Q0E-Q1  UNITS  INCH. 
V-3CRLF.----2.00F  +  00  UNITS  INCH. 


Figure  (5-1) 


Plot  of  total  deviation  of  second, 
third  and  sixth  order  notch-filter 
between  a=0.5  and  a=0.999. 
The  region  in  circle  is  exaggerated 
in  Fig.  (5-2). 
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unimodal  objective  function  of  single  parameter.   The  method 
used  here- -referred  to  by  some  authors  as  the  uniform  or 
factorial  search- -compares  formatly  with  the  Fibonacci 
Search  in  that  it  takes  more  computation  time  than,  but  in 
this  case  of  error-performance  E(a),  in  the  author's 
opinion,  is  the  most  efficient  one.   In  spite  of  the  cost 
of  computation  time,  the  method  is  simple  and  the  result 
is  precise. 

The  results  carried  out  in  this  chapter  agree  with 
computer  experiments  in  Chapter  IV  by  considering  the 
numerical  value  in  Tables  (4-1),  (4-2),  and  (4-3).   Higher 
orders  of  notch- filter  caused  deeper  notch-gains  and  smoother 
response  in  passband.   It  is  noted  that  the  smoother  pass- 
band  is  correct  for  coefficient  close  to  one.   It  is  seen 
in  Fig.  (5-2)  for  noise-performance  in  the  region  close 
to  one. 

The  value  of  the  optimum  coefficient,  a,  increases  with 
the  order  of  filter  provided  the  conditions  (i)  and  (ii) 
are  fixed.   The  higher  order  of  notch- filter ,  the  closer  a 
approaches  one,  the  constraint  is  placed  on  a  as  mentioned 
in  section  D  of  Chapter  II,  due  to  the  finite  wordlength, 
e.g.,  the  single  precision  on  IBM  360/70  has  32-bit 
wordlength.   The  inequality  (2-23)  gives 


a2  =  A  >  2"31 


or 


a2  <  1  -  2~31  ~  0.999999999  (5.5) 


119 


•  ♦ 
.3 


Sj 


•A 


u 


I 
I — \ 


Sixth  order 
Third  order 
Second  order 


j:2 


33; 


X-SCflLr--2.0QF.-Qi   LIMITS   INCH, 
V-SC RLf :•-?..  QGE-0 3  UNITS   INCH. 


C-Q 


a 


Figure  (5-2)   Noise  performance  in  the  region 
of  coefficient  close  to  one. 


1-2  0 


-  31 
The  value  (1-2    )  is  the  upper  limit  of  coefficient,  a, 

in  theory.   In  practical  view,  the  upper  limit  is  smaller 

-  31 
than  (1-2    )  by  the  reason  of  error  introduced  in 

arithmetic  operation  and  error  in  representation  as  stated 

above.   Besides,  in  the  implementation  of  notch-filter, 

the  optimum  coefficient  corresponding  the  larger  wordlength 

in  register  of  computer  will  be  closer  to  unity  than 

the  shorter  one.   For  example,  in  implementation  of  sixth 

order  of  single  precision  on  IBM  360,  the  upper  limit  of 

coefficient  is  0.999,  for  coefficient  exceeds  this  limit, 

the  system  begins  to  degenerate.   It  is  the  same  limit  for 

second  and  third  order.   One  can  see  this  phenomenon 

numerically  in  Tables  (4-1),  (4-2),  and  (4-3)  and  the 

behavior  of  changing  in  steady-state  frequency  response 

plots  of  second,  third  and  sixth  order  in  previous  chapter. 

One  sees  that  for  a  =  0.9999,  the  notch-gain  less 

than  3  dB  in  magnitude  and  the  passband  ripple  is  unacceptable 

It  is  foolish  to  accept  them  as  the  notch- filters . 

31 
Furthermore,  for  a  >  (1  -  2   )  the  system  may  become 

the  oscillator  due  to  error  in  rounding  since  it  could 

make  the  poles  of  their  transfer  functions  exactly  on  the 

unit  circle  of  z-plane. 
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VI.   CONCLUSION 

The  filter  design  was  used  to  design  the  digital  notch- 
filters.   The  effect  of  order  of  filter  on  notch-gain, 
notch-width,  and  passband  ripple  was  investigated. 

From  results  of  this  study,  the  simplest  configuration 
of  digital  notch-filters  was  found;  hence  implementation 
will  be  economical  in  hardware.   The  notch- frequency ,  the 
notch-gain  and  the  notch-width  are  very  easily  controllable 
by  varying  only  one  coefficient  in  realization  as  well  as 
the  agreement  between  the  predicted  and  practice  sampling 
time.   This  property  is  usable  in  an  array  of  notch-filters 
to  eliminate  interference  signals  of  variety  of  frequencies. 

It  was  seen  that  with  selected  coefficient  the  pass- 
band  is  very  smooth  (for  sixth  order  less  than  0.05  dB) 
and  the  notch-width  is  very  narrow  (less  than  one  hertz) . 
The  improvement  can  be  achieved  by  increasing  the  order  of 
filters  provided  the  coefficient  close  enough  to  the  optimum 
coefficient  which  was  found  by  Search  procedure.   The  opti- 
mum computer  program  was  used  and  the  result  was  precise. 
This  program  required  large  computer  time.   That  is  the 
reason  some  authors  do  not  like  it. 
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